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The firing squad synchronization problem (FSSP) on cellular automata has been studied extensively 
for more than forty years, and a rich variety of synchronization algorithms have been proposed for 
not only one-dimensional arrays but two-dimensional arrays. In the present paper, we propose a 
simple recursive-halving based optimum-time synchronization algorithm that can synchronize any 
rectangle arrays of size m x n with a general at one corner in m+n + max(m.n) — 3 steps. The 
algorithm is a natural expansion of the well-known FSSP algorithms proposed by Balzer [1967], 
Gerken [1987], and Waksman [1966] and it can be easily expanded to three-dimensional arrays, even 
to multi-dimensional arrays with a general at any position of the array. The algorithm proposed is 
isotropic concerning the side-lengths of multi-dimensional arrays and its algorithmic correctness is 
transparent and easily verified. 

1 Introduction 

We study a synchronization problem that gives a finite-state protocol for synchronizing large-scale cellu- 
lar automata. The synchronization in cellular automata has been known as a firing squad synchronization 
problem (FSSP) since its development, in which it was originally proposed by J. Myhill in Moore [1964] 
to synchronize all parts of self-reproducing cellular automata. The problem has been studied extensively 
for more than forty years [1-23], and a rich variety of synchronization algorithms have been proposed 
for not only one-dimensional (ID) arrays but two-dimensional (2D) arrays. The ID FSSP is described as 
follows: given a one-dimensional array of n identical cellular automata, including a general at one end 
that is activated at time t = 0, we want to design the automata such that, at some future time, all the cells 
will simultaneously and, for the first time, enter a special firing state. 
Some questions may arise: 

- How can we synchronize multi-dimensional arrays? 

- Can we expand those 2D FSSP algorithms proposed so far to 3D arrays, or more generally to 
multi-dimensional arrays? 

- Can we generalize those 2D FSSP algorithms to generalized ones, in which an initial general is 
located at any position of the array? 

- What is a lower bound of time steps needed for synchronizing multi-dimensional arrays with a 
general at one corner? 

- What is a lower bound for the generalized case? 
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In the present paper, we attempt to answer these questions by proposing a new, simple recursive- 
halving based synchronization algorithm for 2D rectangle cellular automata. The algorithm can synchro- 
nize any 2D rectangle array of size m x n with a general at one corner in m + n + max(m, n) — 3 steps. 
An implementation in terms of local transition rules is also given on a 2D cellular automaton, not only 
for the array with a general at one corner but a generalized case. 

The algorithms proposed in this paper are interesting in the following view points. 

- The 2D algorithm proposed is isotropic with respect to shape of a given rectangle array, i.e. no 
need to control the FSSP algorithm for longer-than-wide and wider-than-long input rectangles. 

- The algorithm proposed can be easily expanded to 3D arrays, even to multi-dimensional arrays. 

- The algorithm is a natural generalization of the well-known ID FSSP algorithms developed by 
Waksman [1966], Balzer [1967] and Gerken [1987]. It gives us a new view point of those classical 
ID FSSP algorithms based on recursive-halving. 

- The algorithm can be expanded to a generalized FSSP solution, where an initial general is at an 
arbitrary position of a given array. 

- The algorithm can be generalized to an optimum-time generalized FSSP solution. 

In Section 2 we give a description of the 2D FSSP and review some basic results on 2D FSSP algo- 
rithms. Section 3 defines the recursive-halving marking on ID arrays and gives some preliminary lemmas 
for the construction of 2D FSSP algorithms. In Sections 4, 5, and 6 we present a new 2D FSSP algorithm 
based on the recursive-halving marking and several multi-dimensional expansions. Two implementations 
in terms of 2D cellular automata are also presented for the optimum-time FSSP algorithms. Most of the 
descriptions of the multi-dimensional FSSP algorithms are based on the 2D FSSP algorithms. Some 
expanded and generalized theorems for multi-dimensional arrays are given without proofs. 

2 Firing Squad Synchronization Problem on Two-Dimensional Arrays 

Figure 1 shows a finite two-dimensional (2D) cellular automaton consisting of m x n cells. Each cell 
is an identical (except the border cells) finite-state automaton. The array operates in lock-step mode in 
such a way that the next state of each cell (except border cells) is determined by both its own present 
state and the present states of its north, south, east and west neighbors. All cells (soldiers), except the 
north-west corner cell (general), are initially in the quiescent state at time t = with the property that 
the next state of a quiescent cell with quiescent neighbors is the quiescent state again. At time t = 0, the 
north-west corner cell Cn is in the fire-when-ready state, which is the initiation signal for the array. The 
firing squad synchronization problem is to determine a description (state set and next-state function) for 
cells that ensures all cells enter the fire state at exactly the same time and for the first time. The tricky 
part of the problem is that the same kind of soldier having a fixed number of states must be synchronized, 
regardless of the size m x n of the array. The set of states and next state function must be independent of 
m and n. 

The problem was first solved by J. McCarthy and M. Minsky who presented a 3«-step algorithm 
for ID cellular array of length n. In 1962, the first optimum-time, i.e. (2n — 2)-step, synchronization 
algorithm was presented by Goto [1962], with each cell having several thousands of states. Waksman 
[1966] presented a 16-state optimum-time synchronization algorithm. Afterward, Balzer [1967] and 
Gerken [1987] developed an eight-state algorithm and a seven-state synchronization algorithm, respec- 
tively, thus decreasing the number of states required for the synchronization. Mazoyer [1987] developed 
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Figure 1 : A two-dimensional (2D) cellular automaton. 



a six-state synchronization algorithm which, at present, is the algorithm having the fewest states for ID 
arrays. 

On the other hand, several synchronization algorithms on 2D arrays have been proposed by Beyer 
[1969], Grasselli [1975], Shinahr [1974], Szwerinski [1982], Schmid [2003], Schmid andWorsch [2004], 
Umeo, Maeda, Hisaoka and Teraoka [2006], and Umeo and Uchino [2008]. It has been shown indepen- 
dently by Beyer [1969] and Shinahr [1974] that there exists no 2D cellular automaton that can synchro- 
nize any 2D array of size m x n in less than m+n + max(m,n) — 3 steps. In addition they first proposed 
an optimum-time synchronization algorithm that can synchronize any 2D array of size m x n in opti- 
mum m + n + max(m,7i) — 3 steps. Shinahr [1974] gave a 28-state implementation. Umeo, Hisaoka 
and Akiguchi [2005] presented a new 12-state synchronization algorithm operating in optimum-step, 
realizing a smallest solution to the rectangle synchronization problem at present. 

As for the time optimality of the 2D FSSP algorithms, the following theorems have been shown. 

Theorem l Be y er I 1969 L Shinahr [1974] fhej-g exists no cellular automaton that can synchronize any 2D array 
of size rnxjiin less than /;; • // • max(m, n) — 3 steps, where the general is located at one corner of the 
array. 

Theorem 2 Shinahl I 1974 !- Umeo < Hisaoka ' and Aki § uchi I 2005 ! There exists a cellular automaton that can synchro- 
nize any 2D array of size mxn at exactly m + n + max(m, n) — 3 steps, where the general is located at 
one corner of the array. 



3 Recursive-Halving Marking 

In this section, we develop a marking schema for ID arrays referred to as recursive-halving marking. 
The marking schema prints a special mark on cells in a cellular space defined by the recursive-halving 
marking. The marking itself is based on a ID FSSP synchronization algorithm. It will be effectively 
used for constructing multi-dimensional FSSP algorithms operating in optimum-time. 

Let S be a ID cellular space consisting of cells Q, Q+i, Cj, denoted by where j > i. Let 

\S\ denote the number of cells in S, that is \S\ = j — i + 1. A center cell(s) C x of S is defined by 
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x= \(i+j)/2 odd 
\(i + j-l)/2,(i + j + l)/2 |5|:even. 

The recursive-halving marking for a given cellular space S = [l...n] is defined as follows: 

Recursive-Halving Marking: RHM 



Algorithm RHM(S) 
begin 

if \S\ > 2 then 

if \S\ is odd then 

mark a center cell C x in 5; 
S L := [l...x]; Ss:= [x...n]; 
RHM L (5 L ); RHM R (S ff ); 
else 

mark center cells C v and C x +\ in S; 
S L :=[l...x];S R :=[x+l...n]; 
RHM L (S L ); RHM R (5 ff ); 

end 



Left-Side Recursive-Halving Marking: RHMl 

Algorithm RHM L (S) 
begin 

while \S\ > 2 do 
if | S\ is odd then 

mark a center cell C x in S; 
S L := [l..jt]; RHM L (5 L ); 
else 

mark center cells C v and C r +i in S; 
S L := [l..jc]; RHM L (5i); 



Right-Side Recursive-Halving Marking: RHMr _ 

Algorithm RHM R (5) 
begin 

while |5| > 2 do 
if | S\ is odd then 

mark a center cell C x in S; 
S R := [x...n]; RHM R (S R ); 
else 

mark center cells C x and C x +\ in 5; 
5 fi := [x +!...»]; RHM R (5«); 



For example, we consider a cellular space S = [1...15] consisting of 15 cells. The first center cell is 
Cg, then the second one is C4, C5 and Cn, C12, and the last one is C2, C3, C13, C14, respectively. In case 
S = [1 ... 17] , we get C9, C5, C13, C3, C15, and C2, Ci6 after four iterations. 

Figure [2] (left) shows a space-time diagram for the marking. At time t = 0, the leftmost cell Ci gener- 
ates an infinite set of signals wi,W2, .., each propagating in the right direction at 1/(2* — 1) speed, 
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Figure 2: Space-time diagram for recursive-halving marking on ID array of length n (left) and some 
snapshots for the marking on 42 (middle) and 71 (right) cells, respectively. 



where k = 1,2,3,..., . The 1/1-speed signal w\ arrives at C n at time t = n— 1. Then, the rightmost cell 
C„ also emits an infinite set of signals W\,W2, .., each propagating in the left direction at 1/(2* — 1) 

speed, where k = 1,2,3,..., . The readers can find that each crossing of two signals, shown in Fig. 2 
(left), enables the marking at middle points defined by the recursive-halving. A finite state realization for 
generating the infinite set of signals above is a well-known technique employed in Balzer [1967], Gerken 
[1987], and Waksman [1966] for the implementations of the optimum-time synchronization algorithms 
on ID arrays. 

We have developed a simple implementation of the recursive-halving marking on a 13-state, 314-rule 
cellular automaton. In Fig. 2 (middle and right) we present several snapshots for the marking on 42 and 
71 cells, respectively. Thus we have: 

Lemma 3 There exists a ID 13-state, 314-rule cellular automaton that can print the recursive-halving 
marking in any cellular space of length n in In — 2 steps. 

An optimum-time complexity In — 2 needed for synchronizing cellular space of length n in the classi- 
cal WBG-type (Waksman [1966], Balzer [1967], and Gerken [1987]) FSSP algorithms can be interpreted 
as follows: Let S be a cellular space of length n = 2n\ + 1, where n\ > 1. The first center mark in 5 is 
printed on cell C„ 1+ i at time ?io-center = 3ni. Additional n\ steps are required for the markings thereafter, 
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yielding a final synchronization at time fio-opt = 3ni +n\= 4n\ =2n — 2. In the case n = 2n\, where 
n\ > 1, the first center mark is printed simultaneously on cells C m and C„ 1+ i at time ?iD-center = 3«i — 1. 
Additional n\ — 1 steps are required for the marking and synchronization thereafter, yielding the final 
synchronization at time fio-opt = 3ni — 1 +n\ — 1 = An\ — 2 = 2n — 2. 

J 3m |5|=2m + l, 

?lD-center = < (2) 

I 3«i — 1 |o| = 2?2i- 

Thus, additional fiD-sync steps are required for the synchronization for a cellular space with the 
recursive-halving marks: 

\n x |5|=2bi + 1, 

?1D- S ync = \ . . (3) 

\n\ — 1 \S \ = 2n\. 




Figure 3: Space-time diagram for synchronizing a cellular space with recursive-halving marking (left) 
and some snapshots for the synchronization on 17 (middle) and 32 (right) cells, respectively. 



In this way, it can be easily seen that any cellular space of length n with the recursive-halving marking 
initially with a general on a center cell or two generals on adjacent center cells can be synchronized in 
\n/2~\ — 1 optimum-steps. In Fig. |5} we illustrate a space-time diagram for synchronizing a cellular space 
with recursive-halving marking (left) and some snapshots for the synchronization on 17 (middle) and 32 
(right) cells, respectively. Thus we have: 

Lemma 4 Any ID cellular space S of length n with the recursive-halving marking initially with a gen- 
erals) on a center cell(s) in S can be synchronized in \n/2~\ — 1 optimum-steps. 

As was seen, the first marking of center cell(s) plays an important role. We print a special mark for 
the first center cell(s) of a given cellular space. On the other hand, for the center cells generated thereafter 
are marked with a different symbol from the first one. 





StepO 



Step 1 



Step 2 



Step 3 



Figure 4: Synchronization schema for 2D cellular automaton. 
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4 An Optimum-Time 2D FSSP Algorithm st x 

4.1 Overview of the Algorithm g$\ 

We assume that an initial general G is on the north-west corner cell Cn of a given array of size m x 
n. The algorithm consists of three phases: a marking phase, a pre-synchronization phase and a final 
synchronization phase. An overview of the 2D synchronization algorithm s/\ is as follows: 

Step 1. Start the recursive-halving marking for cells on each row and column, find a center cell(s) of 
the given array, and generate a new general(s) on the center cell(s). Note that a crossing(s) of the center 
column(s) with the center row(s) is a center cell(s) of the array. 

Step 2. Pre-synchronize the center column(s) using Lemma 4, which is initiated by the general in Step 
1. Every cell on the center column(s) acts as a general at the following Step 3. 

Step 3. Synchronize each row using Lemma 4, initiated by the general generated in Step 2. This yields 
the final synchronization of the array. 




Figure 5: Space-time diagrams for the synchronization algorithm on the 1st, ith, and mth rows of a 
longer-than-wide rectangle array of size mxn, respectively. 



Figure |4] illustrates the synchronization schema for 2D cellular automaton. We assume that m = 
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2m\ + l,n = 1n\ + 1, where mi,n\ > 1. The algorithm operates as follows: 

1. At time t = an initial general on the north-west corner emits a 1/1-speed signal along the first 
row and column to print recursive-halving marks. Once a center mark is printed, it is copied to the 
adjacent row and column. At time t = 3m\ + n\ , a center mark of the center column of the array is 
marked, and the center mark of the center row is marked at time t = 3«i + mi. The center of the 
array is marked at time t = ?2D-center = max(3mi +«i,3ni +mi). 

2. Using Lemma 4, the center column will be synchronized with a tentative pre-firing state at time 

t = ?2D-center+'«l- 

3. Once the center column could be synchronized with the pre-firing state, then the cell C, if , 1+ i ini- 
tiates the synchronization for the ith row for each i such that 1 < i < m. Using Lemma 4, for 
any i, 1 < i < m, the ith row will be synchronized at time t = ?2D-center + mi + «i = max(3mi + 
«i,3ni +mi) +m\ +n\ = max(2rai + \,2n\ + 1) + 2mi +2n\ — 1 = m + n + max(m,?i) — 3. Thus, 
the array can be synchronized at time t = m + n + max(m, n) — 3 in optimum-steps. 
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t = 2u-2 
t - m+n+nt-2 
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Figure 6: Space-diagrams of the synchronization algorithm on the 1st, ith, and mth rows of a wider- 
than-long rectangle array of size mxn, respectively. 

Note that the signal propagation for the recursive-halving marking and the wake-up signal for the 
synchronization are made by the same 1/1 speed. Thus, the synchronization can be performed suc- 
cessfully for each column and the array can be synchronized in optimum steps. Figures [5] and [6] illus- 
trate a space-time diagram for the recursive-halving marking and the synchronization operations on the 
1st, ith, and mth row of a longer-than-wide and wider-than-long array, respectively. One can see that 
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each marking operation has been finished before the arrival of the first wake-up signal for the synchro- 
nization. The algorithm operates in optimum-steps in a similar way for the rectangles such as case 1 : 
m = 2m\ + l,n = 2n\, case 2: m = 1m\,n = 1n\ + 1, and case 3: m = 2m\,n = 2n\. 
Thus, we can establish the following theorem. 

Theorem 5 The synchronization algorithm s£\ can synchronize any m x n rectangular array in optimum 
ra + 7i + max(m,?i) — 3 steps. 

We have implemented the algorithm srf\ on a 2D cellular automaton having 60 states and 13633 local 
rules. In Figures [7] and [8] we present some snapshots of the synchronization processes of the algorithm 
on 9 x 12 and 12x9 arrays, respectively. 



5 Expansion to Multi-Dimensional Arrays 
5.1 Three-Dimensional Arrays 

In this section, we show that there exists no algorithm that can synchronize any 3D array of size mxnxi 
with a general at an arbitrary corner in less than m + n + l + max(m,n,£) — 4 steps. 

Theorem 6 The minimum time in which the firing squad synchronization could occur is no earlier than 
m + n + £ + max(m,n,£) — 4 for any 3D array of size mxnxi with a general at an arbitrary corner cell. 

Proof. The proof is made by contradiction. Without loss of generality, we assume that I <m <n. 
It is assumed that there is a cellular automaton ^ that can synchronize an array of size m^xn^x £ 
(£q < mo < no) at step t = to such that: 
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t < m + 2no+4-4 (4) 

Now consider the state of cell C m „ i £ at time t = to. Let i and k be any integer such that 1 < i < mo 
and 1 < k < £q. Consider the signal propagation from the cell Ci i i to C mo i e via any cell Q „ It 
takes: 



(i-l) + (no-l) + (*-l) + (mo-0 + (no-l) + (4>-*) 

= /n + 2n + 4-4 (5) 

steps for the signal to travel from Ci i i to C,„ i e via any cell C, „ The state of the cell C,„ i e at 
step i = to entered the final firing state unaffected by any cells on the plane {Q < i < mo, 1 < it < lo 
}. Therefore, if another three-dimensional array of size mo x tiq x £ was added to the right side of the 
original array (that is, the new array is of size mo x 2no x Iq), the cell C m i £ would still enter the final 
firing state at step t = to- This is because the cell structure ^# is fixed, cell operation is deterministic 
and nothing has changed as far as the cell C m \ i Q is concerned. Since to < mo + 2no + £q — 4, the cell 
Cm i to wm stm b e i n a quiescent state at time t = to. Therefore the cell structure does not represent 
a solution and this is a contradiction. In a similar way, the argument carries over in the cases such as 
£ <n <m,n < I <m, and so forth. 



□ 
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The synchronization algorithm si\ for 2D arrays can be easily expanded to 3D arrays. See Figure [9] 
which illustrates the synchronization schema for 3D cellular automaton. Thus, we have: 

Theorem 7 There exists an optimum-time synchronization algorithm sfi that can synchronize any three- 
dimensional array of size «|X« 2 x 723 with a general at Ci, 1,1 in optimum n\ + JI2 + 113 +max.(ni,n2,n^) — 
4 steps. 




Step Slep 1 step 2 Step 3 Step 4 



Figure 9: Synchronization scheme for a three-dimensional cellular automaton. 



5.2 Multi-Dimensional Arrays 

A kD FSSP algorithm is sketched as follows: 

Step 1. Start the recursive-halving marking on cells along each dimension, find a center cell(s) of the 
given array, generate a general(s) on the center cell(s), and pre-synchronize the center point(s): zero- 
dimensional sub-array of the array. 

Step 2. Pre-synchronize a ID sub-array along the 1st dimension containing the pre-synchronized center 
cell. 

Step 3. - Step k. For j = 2 to k — 1, by increasing the number of dimensions, pre-synchronize a 7D 

sub-array containing the pre-synchronized (j — 1)D sub-array. 

Step k + 1. Synchronize the kD array. This yields the final synchronization of the given array. 

Theorems 6 and 7 can be expanded to the kD arrays. 

Theorem 9 There exists no cellular automaton that can synchronize any kD array of size n\ x n% x ... x 
with a general at Cij....^ in less than £f =1 n, +max(ni,n2,...,nfe) — k — 1 steps 

Theorem 10 There exists an optimum-time synchronization algorithm that can synchronize any &D 

array of size n\ x 112 x ... x with a general at C^i 1 in optimum J^Lj n,-+tnax(ni,/i2, ...,«#) — k— 1 

steps. 



6 Generalization as to General's Position 
6.1 Generalized FSSP on ID Arrays 

The recursive-halving marking scheme on ID array can be easily expanded to the generalized case where 
the initial general is located at any position of the array. Figure [10] illustrates a space-time diagram for 
the recursive-halving marking on ID array of length n with a general on 1 < k < n. The marking is 
based on the generalized FSSP algorithm proposed by Moore and Langdon [1968]. 
We have seen the following theorems for the ID generalized case. 
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Figure 10: Space-time diagram for recursive-halving marking on ID array of length n with a general at 
any position. 

Theorem 11 Moore and Langdon [1968] m i n j mum time in which the generalized firing squad synchro- 
nization could occur is n — 2 + max(k, n — k + 1) steps, where the general is located on the kth cell from 
left end. 

Theorem 12 Moore and Lan § don I 1968 I There exists a 17-state cellular automaton that can synchronize any 
one-dimensional array of length n in optimum n — 2 + max(k, n — k+l) steps, where the general is located 
on the kth cell from left end. 

An optimum-time complexity n — 2 + max(k, n — k+l) needed for synchronizing cellular space of 
length n in the classical ML-type (Moore and Langdon [1968]) generalized FSSP algorithm can be 
interpreted as follows: Let S be a cellular space of length n = 2n\ + 1, where n\ > 1. The first center 
mark in S is printed on cell C„ 1+ i at time fmg-center = 3«i — min(& — l,n — k). Additional n\ steps are 
required for the markings thereafter, yielding a final synchronization at time ?iDg-opt = 3«i —mm{k — 
\,n — k)+n\ = \n\ — min(& — I ,n — k) = 2n — 2 — min(& —l,n — k) = n — 2 + max(k,n — k+ 1). 

In the case n = 2n\, where n\ > 1, the first center mark is printed simultaneously on cells C m and 
Qu+i at time fmg-center = 3«i — 1 — min(& — l,n — k). Note that two cells C m and C„ 1+ i are pre- 
synchronized at time ?iDg-center = 3«i — 1 — min(fc — l,n — k). Additional n\ — 1 steps are required for 
the marking thereafter, yielding the final synchronization at time £iDg-opt = 3«i — 1 — min(& — 1 , n — k) + 
72 1 — 1 = 4rci —2 — min(& — I ,n — k) = 2n — 2 — min(/c — l,n — fc) = n — 2 + ma.x(k,n — k + 1). 



j 3«i — min(& — l,n — fc) |£|=2ni + l, 

flDg-center = < , . (6) 

I 3«i — 1 — mm(£ — — k) \S\=2n\. 
In addition, fiDg-sync steps are required for the synchronization for a cellular space with the recursive- 
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halving marks: 

J«i \S\=2ni + l, 

flDg-sync - < . . ( /) 

I «i — 1 |o| = 2ni. 

6.2 Generalized FSSP on 2D Arrays 

By a similar method employed in Section 4, we can develop the following theorem for the generalized 
case. 

Theorem 13 There exists no 2D cellular automaton that can synchronize any 2D array of size mxn with 
an initial general on C r , in less than m + n + max(m,n) — min(r,ra — r+ 1) — mm{s,n — s+ 1) — 1 steps, 
where 1 < r < m, 1 < s < n. 

Now we are going to present a generalized optimum-time FSSP Algorithm 3/4 for 2D arrays. We 
assume that an initial general G is on the cell C r}S of a given array of size mxn, where 1 < r < m, 1 < 
s < n. The algorithm consists of three phases: a marking phase, a pre-synchronization phase and a final 
synchronization phase. An overview of the 2-D synchronization algorithm ^ is as follows: 

Step 1. Start the recursive-halving marking for cells on each row and column, find a center cell(s) of 
the given array, and generate a new general(s) on the center cell(s). Note that a crossing(s) of the center 
column(s) with the center row(s) is a center cell(s) of the array. 

Step 2. Pre-synchronize the center column(s) using Lemma 4, which is initiated by the general in step 
1. Every cell on the center column(s) acts as a general at the next Step 3. 

Step 3. Synchronize each row using Lemma 4, initiated by the general generated in Step 2. This yields 
the final synchronization of the array. 

The array can be synchronized at time t = m+n+ranx{m,n) — min(r,m — r+ 1) — mm{s,n — s+ 1) — 1 
in optimum-steps. 

Theorem 14 There exists an optimum-time synchronization algorithm ^ that can synchronize any mxn 
rectangular array with a general at C r . v in optimum m + n + max(m, n) — min(r, m — r + 1) — min(^, n — 
s + 1) — 1 steps, where l<r<m, I <s <n. 

We have implemented the algorithm ^ on a 2D cellular automaton having 269 states and 163662 
local rules. We have checked the rule set for any array of size mxn, with 2 < m,n < 100, and any 



general's position in the array. In Figure 1 1 we present some snapshots of the synchronization processes 



of the algorithm on an 8 x 13 array with a general on C3 5. 



6.3 Generalized FSSP on Multi-Dimensional Arrays 

Theorems 13 and 14 can be expanded to three or more dimensional arrays. 

Theorem 15 The minimum time in which the firing squad synchronization could occur is no earlier than 
n\ +722 +«3 + max(«i ,«2,«3) — min(ri,7ii — r\ + 1) — min(r2,7i2 — ?"2 + 1) — min(r3,«3 — r^ + 1) — 1 for 
any 3D array of size fi[X« 2 x «3 with a general at C run ^. 

Theorem 16 There exists an optimum-time synchronization algorithm .2^5 that can synchronize any 
3D array of size n\ x « 2 x «3 with a general at C ri .,- 2jr3 in optimum n\ +«2 + «3 + max(ni,n2,«3) — 
min(ri ,«i — n + 1) — min(r2,?i2 — r2 + 1) — min(r3,?i3 — r^ + 1) — 1 steps. 
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Theorem 17 There exists no cellular automaton that can synchronize any kD array of size n\ x ?i2 x ... x 
nk with a general at C n ,r 2 ,...,r k in l ess than Yli=\ n/+max(ni, 712, — L?=i rniri(r;, n,- — r,- + l) — 1 steps. 

Theorem 18 There exists an optimum-time synchronization algorithm jz^ that can synchronize any kD 

array of size n\ x ni x ... x with a general at C ri , r2 n in optimum Lf =1 «; + max(ni,n2, ...,nfe) — 

£* =1 min(r,,ra; — r,- + 1) — 1 steps. 



7 Conclusions 

We have proposed a new class of optimum-time multi-dimensional FSSP algorithms based on recursive- 
halving marking. The class includes the well-known optimum-time FSSP algorithms developed by 
Waksman [1964], Balzer [1966] and Gerken [1987] with a general at one end and Moore and Lang- 
don [1968] with a general at any position. 
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